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ABSTRACT 



We report on the results of an optical spectroscopic survey designed to con- 
firm the youth and determine the spectral types among a sample of young stellar 
object (YSO) candidates in the Serpens Molecular Cloud. We observed 150 
infrared excess objects, previously discovered by the Spitzer Legacy Program 
"From Molecular Cores to Planet- Forming Disks" (c2d), bright enough for sub- 
sequent Spitzer /IRS spectroscopy. We obtained 78 optical spectra of sufficient 
S/N for analysis. Extinctions, effective temperatures and luminosities are esti- 
mated for this sample, and used to construct H-R diagrams for the population. 
We identified 20 background giants contaminating the sample, based on their 
relatively high extinction, position in the H-R diagram, the lack of Ha emission 
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and relatively low infrared excess. Such strong background contamination (25%) 
is consistent with the location of Serpens being close to the Galactic plane (5° 
Galactic latitude). The remaining 58 stars (75%) were all confirmed to be young, 
mostly K and M-type stars that are presumed to belong to the cloud. Individ- 
ual ages and masses for the YSOs are inferred based on theoretical evolutionary 
models. The models indicate a spread in stellar ages from 1 to 15 Myr, peaking 
at 2 - 6 Myr, and a mass distribution of 0.2 to 1.2 Mq with median value around 
0.8 Mq. Strong Ha emission lines (EW[Ha] > 3 A) have been detected in more 
than half of the sample (35 stars). The mass accretion rates as derived from 
the Ha line widths span a broad distribution over 4 orders of magnitude with 
median accretion rate of 10~^ Mq yr~^. Our analysis shows that the majority 
of the infrared excess objects detected in Serpens are actively accreting, young 
T-Tauri stars. 

Subject headings: ISM: individual (Serpens) - stars: pre-main sequence - stars: 
Hertzsprung-Russell diagram - stars: planetary systems: protoplanetary disks 



Introduction 



The Serpens molecular cloud has received considerable attention over the past decade. Be- 
cause it is an actively star-forming complex containing a substantial mass of molecular 
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The Serpens cloud is one of the five clouds selected as part of the Spitz er Legacy program 
"From Molecular Cores to Planet-Forming Disks" (c2d; lEvans et al.ll2003l ). providing images 
in the 3-70 /im range. The wide wavelength coverage and high sensitivity in the infrared 
(IR) of Spitzer make it possible to easily identify a new, complete, flux limited (down to 
luminosities below 0.01 Lq) young stellar population. This population has been found to be 
distributed over almost the entire area surveyed, offering an opportunity to determine the 
stellar content in different regions of the cloud, the distributions of the youngest stars and 
substellar objects, and the properties of their circumstellar envelopes and disks. 

The c2d prog ram has mapped a 0.89 deg^ portion out of the more than 10 deg^ area 



(Kaas et al 



2004r) o f the Serpens molecular cloud. Assuming a distance of 259 ± 37 pc 



( IStraizys et al.lll996l ). this corresponds to a covered area of about 2.5 x 9 pc. This region 
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was discovered to be very rich in young stars (jHarvey et al.ll2007bl ). some of them previously 
identified with ISO data (iDjupvik et al.ll2006l ). This poorly known cluster of young stars, 
cluster B, and the previously unknown star-forming region around it, is located about half 



a deg ree southwest of the well-known Serpens Cloud Core, containing cluster A (IKaas et al. 



20041 ) ■ It has a high density of young stars, making it an unique target region for obtaining a 
full and consistent picture of clustered low-mass star formation and compare this with young 
stars in the surrounding field. 

We are currently carrying out a mult i- wavelength survey of this region, from X-ray to 
millimeter wavelengths, in order to create a "template" sample for the study of the evolu- 
tion of circumstellar disks around stars younger than ~ 10 Myr, within a single, small and 
well defined region. This work provides the necessary information for anchoring the study of 
protoplanetary disks to their parent population by means of the optical spectroscopic classifi- 
cation of the central stars. Precise stellar properties are needed to perform an accurate study 
of the evolutionary stages of the young stellar objects and their disks. Most spectroscopic 
studies of protoplanetary disk evolution refer to samples of young stars scattered across the 
sky or to sources in large star-forming clouds like Taurus, making it difficult to separate 
intrinsic evolutionary effects from those caused by external influences such as environment 
or star formation history. 

In this paper we report on our optical spectroscopic survey designed to conflrm the youth 
and determine spectral types of the newly discovered young stellar object (YSO) candidates 
in the Serpens region with Spitzer. These spectra will also be used to es t imate mass accretion 



rates using the st rength of the Ha emission line ( iMuzeroUe et al.ll2003l : IWhite fc Basril 12003 
Natta et al.ll2004j ). and the extinction toward each object. 



The paper is constructed as follows: § [2] describes the selection criteria for our sample, 
and in § [3] we describe the observations and data reduction. In § H] we discuss our spectral 
classiflcation methods, along with its results for this sample and derive effective temperatures. 
In § O we present our extinction studies based on the Spitzer photometry and on our optical 
spectra. In §[6] we derive luminosities for the objects. This allows us to place each object in 
an Hertzsprung-Russell (HR) diagram, and derive its age and mass from comparison with 
evolutionary models. The Ha line and accretion rates are discussed in § [71 Finally, in § [H] 
we present our conclusions. 
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2. Sample Selection 



From the c2d catalog, iHarvey et al.l (l2007bl ) plotted spectral energy distributions (assuming 
a K7 spectral type for all the sources) and identified 235 YSO candidates - objects with 
IR excess, interpreted as due to disks or envelopes. These authors discuss the problem that 
background galaxies and bright background post-AGB stars may show the same color distri- 
bution as that of young stellar objects. Particularly, the post-AGB stars could contaminate 
the bright end of the luminosity distribution. 

To this original sample of 235 objects, a fiux threshold of 3mJy at 8fim was applied, 
resulting in 150 sources bright enough to be observed in follow-up observations with the 
InfraRed Spectrograph (IRS), onboard Spitzer. The IRS observations, carried out to study 
the evolution of dust in pr otoplanetary disks of a complete, fiux limited sample will be 



presented in a future paper (lOliveira et al.llin prep.l ). Here, we present optical spectroscopy 



of 78 stars bright enough to obtain optical spectroscopy with a 4m-class telescope. The 
remaining objects are too extincted in the optical to be detected, and are being classified in 
an ongoing near-IR spectroscopic survey. 

The resulting sample contains objects which span a range in evolutionary stages from 
embedded Class I (i.e., spectral slope Q;2-24^jm > 0.3) sources to young stars with protoplan- 
etary disks (Class II sources, -1.6 < a2-24/im < 0.3) and stars with photosph eric fiuxes in 



the I RAC bands and excess in the MIPS bands (also called "cold disks", e.g. iBrown et al. 



20071 ) 



Figure [T] shows the spatial distribution of the complete sample of 150 Serpens YSOs ob- 
served with Spitzer/IRS. The filled and open symbols correspond to detected and undetected 
objects from the optical survey presented in this paper, respectively. 



3. Observations and Data Reduction 

The spectroscopic data were obtained during observing runs at three telescopes: two runs 
with the multi-object Wide Field Fibre Optical Spectrograph {WYFFOS) at the 4.2m 
William Herschel Telescope ( WHT) in 2006 May and 2007 June, and one run with the 
3.6m Device Optimized for the LOw RESolution (DOLORES) at the Telescopio Nazionale 
Galileo (TNG) in 2006 June, both telescopes are located at the Observatorio del Roque de 
los Muchachos in La Palma, and one run with the Calar Alto Faint Object Spectrograph 
(CAFOS) at the 2.2m Calar Alto Telescope. No absolute fiux calibration was performed, nor 
were telluric absorption features removed. Technical information on each run is specified in 
Tabled! Figure [2] shows a representative sample of the spectra of the objects in our sample. 
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Fig. 1. — Observed objects shown on the 8.0/im IRAC map of the Serpens Molecular Cloud, 
imaged as par t of the Spitzer Leg acy Program "From Molecular Clouds to Planet- Forming 
Disks" (c2d) (IHarvey et al.ll2006l ). Filled symbols are the 78 objects for which we present 
optical spectroscopy in this paper: filled circles are confirmed members of the cloud and 
filled triangles are likely candidate background contaminants (§[6]). The open symbols are 
the remaining young stellar object candidates in the sample, not detected in our survey of 
150 YSOs. Overlaid are the contours (5, 10, 15, 20 and 25 mag) of visual extinction derived 
from the "c2d" extinction maps. 
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3.1. WHT data 



WYFFOS wa.s fed from the AutoFib 2 robotic fibre positioner (iParry et al.lll994j ) and targets 
were selected to maximize the number of objects observed in each configuration, using its 
program AF2_C0NFIGURE. We used the "small" fibre module of WYFFOS, corresponding to 
fibre aperture of 1.6". The wavelength range was ~3000 A, centered at 7000 A, in order 
to include the Ha and Li I (6707 A) lines. The average resolving power was R ~ 1750. 
In order to adequately remove the effects of cosmic ray hits on the detector, each field was 
observed using three exposures. Flat-fields and bias frames were obtained at the beginning 
of each night and neon calibration arc lamp spectra at the beginning of the nights, as well 
as separately for each configuration. Unused fibres were placed on the sky. We obtained 
WYFFOS spectra for 71 objects of our sample. 

Data reduction was performed using our own code, based on the instrument reduction 
guidelin^, within IRAF and IDL. After bias subtraction and fiat-fielding, the spectra were 
extracted and wavelength calibrated with the arc Ne comparison lamp data, to a precision 
of ~0.5 A. 



3.1.1. Sky subtraction 

The quality of the optical spectra is limited by the accuracy of the sky subtraction, especially 
for faint stars. As opposed to long-slit spectroscopy, where an adequate sky removal can be 
performed using the sky background on opposite sides very close to the object spectrum, a 
precise sky subtraction in multi-object fibre spectroscopy is complicated b y the lack of a sky 



spect rum close to the object spectrum both in position and in time (e.g., see lWyse fc Gilmore 



19921 ). In general, the sky background depends mainly on the local sky conditions at the 
observatory during the time of observations. However, in the case of a molecular cloud, the 
"sky" may also have a significant contribution from extended continuum and line emission 
from the cloud itself. 

We used the fibres placed on "empty" positions in the field to measure the contribution 
from the sky. The sky spectra were reduced and extracted in an identical manner to the 
science data. The sky spectra were visually inspected and a few of bad quality were discarded. 
None of them showed diffuse Ha in emission from the cloud. To reach optimum results, even 
for the faintest sources, a single master sky spectrum, built with an average of 12 - 15 sky 
spectra, was computed for each configuration. During the sky subtraction of each star, we 



http://www.ing.iac.es/Astronomy/instrunicnts/af2/index.html 
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Fig. 2. — Selected spectra of a representative sample of the observed objects. Along with 
their number from Table [Tj the extracted spectral types are shown for reference. The strong 
Ha lines are noticeable in several spectra. The main remaining atmospheric features are 
indicated with ©. 

applied small scalings to the master sky spectrum to take into account the fact that the 
throughput of different fibres shows slight variations. The scalings were typically in the 
range 0.95-1.05, and chosen to optimize the sky subtraction. 

In several cases, the sky subtraction proved imperfect with residuals being most notable 
at A >7700A. These residuals are very common in low to intermediate S/N fibre spectra due 
to the strong OH sky lines beyond 6700 A. It is very difficult to remove these OH features 
because the sky lines are usually not well resolved, and also due to the non-linear variations 
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in throughput and dispersion for different fibres. We have checked, however, that the spectral 
type classification and Ha equivalent width measurements are not substantially hampered 
by these effects. 

3.2. TNG data 

The observations with DOLORES used the long-slit mode. Slit rotation was calculated in 
order to enable more than one object in the slit whenever two targets were closer than the 
length of the slit (8'). The wavelength range of ~1600 A was centered at 7000 A, again 
including the Ha and Li I (6707 A) lines, with resolving power R ~ 2500. Exposure times 
ranged from 12 to 1000 seconds. Flat-fields, bias and argon comparison arc lamp spectra 
were obtained at the beginning of each night. 

After bias subtraction and flat-fielding, the spectra were extracted and wavelength cal- 
ibrated with a combination of the Ar arc lamp data and airglow line^, including the OH 
Meinel bands, using the "doslit" task of the specred package in IRAF. Sky subtraction 
was done during the extraction based on the off-source regions along the slit. In total, we 
obtained DOLORES spectra of 18 sources, 15 of which were also observed with WYFFOS. 

3.3. CAFOS data 

CAFOS was used with mask observations in 2008 June. Three masks were taken using the 
grism R-100. Flat-fields, bias and HgCd/He/Rb comparison arc lamp spectra were obtained 
at the beginning of each night. The wavelength coverage (3500 A, centered on 7750 A) for 
each source varies slightly depending on their position on the CCD, with resolving power 
R ~ 3523. The exposure time of 2400 seconds for each mask was divided in two for cosmic 
ray rejection. 

After bias subtraction and flat-fielding, the fields were divided so that each source would 
be in an image by itself and then further reduced as standard slit observations: the spectra 
were extracted and wavelength calibrated using the "doslit" task of the specred package in 
IRAF, performing sky subtraction during the extraction based on the off-source regions along 
each slit. 

The low resolution of these instruments, combined with low S/N for most faint objects. 



^http: / / astro. u-strasbg.fr/^koppen/divers/AirGlow. html 
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did not allow us to detect or put meaningful upper limits on the Li I (6707 A) line. 



4. Spectral Classification 

Two methods were used in the spectral classification of the sources and are explained 
below. 



4.1. Method I 



The o bjects were classified following the spectral classification scheme of iHernandez et al. 



( I2OO4J ) optimized to classify early type stars (up to early G). For stars with spectral type later 
than G, the scheme was extended by i ncorporating spectra l indices for the ranges G0-K5 
and K5-M6 Jsiciha-Aguilar et all [20051 . ISricefio et aUboodj^h . More information about this 



classification scheme can be found in the references above. In short, it is based on strong 
spectral features that are sensitive to Tcs but insensitive to reddening in a specific luminosity 
class (dwarf in this case), since it uses local continua for the determination of the feature's 
equivalent widths (EW). It is also designed to avoid problems caused by non-photospheric 
contributions (e.g. emission lines). The wavelength range of the WYFFOS spectra covers 12 
of these features for all spectral types, while DOLORES covers 6 and the range of CAFOS 
covers 15 features. On average, 6 features for WYFFOS, 5 for DOLORES and 7 features for 
CAFOS were used for each classification, as exemplified in Figure [31 

The spectral types calculated from each feature result in a weighted mean spectral 
type. The errors in the estimated spectral types have two contributions: the error in the 



measurement of each feature and the error of its fit to the standard (see Hernandez et al. 



2004 for more details). A few objects were observed in both runs (as can be seen in Tabled]). 
For these objects, both spectra were ran separately through the spectral classification code. 
Their spectral classifications were found to be consistent with each other, typically within 
one subclass. 



•^Information about the code at |http:/ /www. astro. lsa.umich.edu/~hernandj/SPTclass/sptclass. html 



- 10 - 




6500 7000 7500 
Wavelength (A) 



Fig. 3. — Features used in the spectral type classification scheme of iHernandez et al.l (120041 ) 
are shown on the spectrum of source #113. 



4.2. Method II 



Due to telluric absorption features not corrected for in the spectra presented here, care needs 
to be taken when using automatic fitting of features very close to these telluric bands. To 
make sure that the spectral types derived in Method I are correct, an additional spectral 
typing procedure was applied to the data. 



The library of standard spectra from Mora et al.l f 20011 ) . covering the wavelength range 
of 5800-7000A, was used in this method. In this region, different lines and bands are found 
for specific spectral types, making them unique. Each science spectrum was normalized and 
compared to normalized standards of different spectral types. The spectral classification 
procedure was based on finding the best correlation between observed and standard spec- 
tra, combined with a visual inspection of each match. Figure H] shows an example of the 
correlation. 



4.3. Results 

The two methods agree in their classification (within 1 subclass) for about 80% of the 
sources. When in disagreement. Method II is judged to be more reliable because of the 
visual inspections, whereas Method I may be contaminated by the forementioned telluric 
absorption features. The spectral types obtained for our sample are given in Table [2l Figure 
|5] shows the distribution of spectral types, which is clearly dominated by K and M-type stars. 
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Fig. 4. — Correlation between science spectrum of source #115 (red) and standard spectra 
for different spectral types (b lack; spectral type s indicated on the top right of each standard 
spectrum from the library of iMora et al.l (l200ll )). 
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Fig. 5. — Distribution of spectral types of YSO candidates in the Serpens Molecular Cloud 
based on the classification scheme described in § HI 
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4.4. Effective Temperature 



Effective temperatures {Tes) were determined using calibratio ns relating tfiem to spectra l 
types. For stars earlier than MO, the relationsh ip established bylKenvon fc HartmannI (119951 ) 
was adopted, while for later type stars that by lLuhman et al.l (120031 ) was used. This method 
was also used in the determination of T ^ff for similar study on another of the clouds observed 
by the c2d, Cha II (jSpezzi et al.ll2008l ). The errors in Teg come directly from the errors in 
the spectral type determination. 



Visual Extinction 



Once the spectral types are kno wn, we used synthetic NEXTG EN spectra (IHauschildt et al. 
19991 ) and the extinction law by IWeingartner fc Draind (120011 ) to compute the observed ex- 
tinction towards our targets (AyOhs). The appropriate synthetic spectrum for each spectral 
type was extincted with Ay ranging from to 20 mag (with increments of 0.5 mag) and a 
minimization was performed between the observed and synthetic extincted spectra to esti- 
mate the extinction. We estimate an intrinsic error in our AyOhs determinations of 1.0-1.5 
mag by compa ring our sou rce 98 (spectral type A3) with the A3 V absolute flux standard 
star Kopff 27 (IStond 119771 ). from the ING list of spectro-photometric standard stars. This 
error includes the uncertainty in the total extinction determination and the fact that the 
fiber-fed spectra are not flux calibrated. 

Table [2] compares our extinction values from the optical spectra (AyObs) with those 
from the "c2d" extinction maps (Ay Cloud). The "c2d" map was constructed using the 
average extinction toward background stars within a 5' beam. Thus, these values give an 
indication of the average amount of extinction that the cloud produces in a given area, 
albeit with a large beam. More information about this extinct ion determination and its 
uncertainties can be found in the "c2d" delivery documentation (lEvans et al.l 120071 ). and in 
Harvev et all J2006h . 



Figure [6] compares the two estimates of extinction for our targets described above. A 
large spread around the line of equal extinction is seen. This dispersion is useful to pinpoint 
objects with peculiar extinction characteristics. Taking into account a typical error of 2 mag 
in both extinction determinations and the large beam sizes for the "c2d" extinction values, 
the larger differences might be explained by two effects: 

(i) For the sources with AyOhs < Ay Cloud, an indication of their depth in the cloud 
can be obtained. Specifically, such sources are likely embedded in the cloud, or located in 
front of it. 
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Fig. 6. — The relationship between the extinction values estimated from our optical spectra 
{AyOhs) and the extinction values from the "c2d" Spitzer data averaged over a 5' radius 
(AyCloud) is shown with crosses. The solid line indicates AyOhs = AyCloud, with typical 
errors of 2 magnitudes (das hed lines). Objects m arked with a square are accreting stars 
according to the definition of I White fc Basril ( 120031 ). while objects marked with a filled circle 
are confirmed post-AGB stars. The shaded region is where we propose background stars 
should lie. 



(ii) If AyOhs » AvCloud (shaded area in Figure [6]), the molecular cloud by itself 
is unlikely to explain such extinction. One one hand, for young stars embedded in their 
parent envelopes, an extra extinction compared to the local average in the cloud is naturally 
expected. Such embedded sources, however, are also expected to show a strong infrared 
excess. On the other hand, stars at distances much larger than that of the cloud are also 
likely to have an enhanced extinction. We argue that this is the dominant category of sources 
in the shaded area of Figure M, though confirmation through another method is needed. 

Indeed, none of the objects in the shaded area was found to have a strong infrared 
excess, as evidenced by the low Spitzer-derived slope (a2-24/im) values, shown in Table [2j 
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This suggests that those objects are not surrounded by large amounts of circumstellar dust 
but rather by tenuous dust layers, consistent with them being post-AGB stars at distances 
of several kpc. Follow-up c2d observations of bright YSO candidates in the Serpens cloud 
with Spitzer IRS confirmed that two of the sources with AyOhs » AyCloud are post-AGB 
stars, seen through the cloud (filled circles in Figure [6]). Post-AGB stars are so bright that, 
at a farther distance, their colors and magnitudes can be comparable to those of YSOs in 
the cloud. They can, however, be distinguished b y the presence of spectral features at 13 
and 21 fj,m, not seen in YSOs (IWaters et al.lll999l ). 



6. H-R Diagram 

Given the spectral type, therefore effective temperature (derived in § 14.41) . and luminosity 
of a source, individual age and mass can be estimated by comparing its position in the 
Hertzsprung-Russell (HR) diagram with evolutionary tracks. 



6.1. Luminosities 

Stellar luminosities were calculated by integrating the NEXTGEN stellar photosphere for 
each object, scaled to dereddened optical (from literature) or 2MASS fluxes, depending on 
their availability (see Figure ?? for two examples). The errors are derived from the uncer- 



ods f 


□r luminosity estimation have been widely used in the literature ( 


van den Ancker et al. 


1997 


van den Ancker et a . 


1998, 


Alcala et al. 


2008. 


Merm et al. 


2008 


). The extinction law 



of IWeingartner fc Draind (120011 ) was used to correct for the reddening, to be consistent with 
the analysis in § [51 A careful inspection of all the resulting SEDs was performed to find 
the best-fit extinction from the two Ay values discussed above plus the best-fit Ay from 
the spectral type and the optical or near-IR photometry. The resulting best-fit extinctions 
(AyFinal) and luminosities are shown in Table [31 Objects # 44 and 47 were not detected by 
2MASS and, therefore, no J, H and K magnitudes are available for them. For this reason, 
we could not determine their luminosities and they are not included in the H-R diagram. 
Objects # 41 and 117 are rising sources (likely due to an edge-on disk or envelope), resulting 
in an unrealistically low luminosity and were marked with a different symbol (square) in the 
H-R diagram. 
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6.2. Results 



The H-R diagram of the sample is s hown in Fi gure [71 Th e stars are compared with evo- 
lutionary tracks from iBaraffe et al.l (Il998l ) and ISiess et al.l (|2000| ). 20 objects are outside 
the range of the isochrones in the HR diagram (gray circles in Figure [7]). These objects are 
too luminous to be at a distance of 259 pc even in the extreme case of having vanishing 
extinction. As discussed above, their SEDs indeed show very little infrared excess, only at 
the longest wavelengths, and therefore they likely do not belong to the cloud. Obviously for 
these objects, no ages or masses could be determined. 

16 of these outliers have AyOhs » AyCloud and therefore appear in the shaded area 
in Figure [61 Thus, there is a 92% correspondence of both methods to identify background 
sources. The 20 background stars are removed from the sample, leaving the remainder as a 
purer YSO sample in Serpens. 

These remaining objects, shown with black circles in Figure [3 are found to be consis- 
tently along the 1 to 10 Myr isochrones and between the 0.2 and 1.2 Mq mass tracks in both 
diagrams, indicating a possibly coeval population of young stella r objects dorn i nated by very 
low-mass stars. Early type (Tpff > 6000K) st ars fall outside the iBaraffe et al.l (Il998l ) tracks, 
but are in the range of the ISiess et al.l (120001 ) models. In these cases, ages and masses could 
only be derived from the latter. 

Individual ages and masses of the YSOs in Serpens wer e derived from the models and 
are shown in Table [31 and Figure [81 The IBaraffe et al.l (Il998[) tracks tend to produce larger 
masses and ages than those found with the Siess et al. (j2000 ) tracks for the same position 
in the diagram. We also find that, for the different models, the differences in ages are within 
the luminosity errors, while the differences in masses are slightly larger than those due to the 
errors. In spite of the differences between models, we find similar mass and age distributions 
in both cases. Both sets of tracks imply a population of YSOs concentrated between 1 and 
15 Myr with a few sources older than 15 Myr, and st rongly peaked at 2 - 6 Myr. The inferred 
median age is 7.5 Myr with the IBaraffe et al.l (119981 ) and 4.7 Myr with the lSiess et al.l (120001 ) 
tracks. Individual m a sses r ange from 0.2 to 2.2 M (^, with median values of 0.8 and 0.6 Mq 
for the IBaraffe et al.l (119981 ) and ISiess et al.l (120001 ) models, respectively. The objects more 
massive than 2.0 Mg (#62, 70, 98, 108, 120, 139, 141) do not lie in the more clustered 
regions of the cloud, and are also separate from each other. The main uncertainty is caused 
by the uncertaiii ty in extinction: higher (lower) extinction values by 2 mag can shift the 
Siess et al.l ( l2000l ) median age by a factor of two to younger (older) ages. It is also important 
to stress that this is not a complete IMF study - we observe a bias towards a larger mean 
age due to the bias in the selection of the sources (i.e. lack of very young, deeply embedded. 
Class objects). Nevertheless, given that some of the stellar ages are clearly larger than 




Fig. 7. — H-R diagram s of the obiects in thi s sample. Overlaid are the iso chrones (blue) and 
mass tracks (red) from lBaraffe et aD ( 1998 ) (left) and Siess et al. ( 2000 ) (right), with their 
corresponding age (in Myr) and mass (in Mq) indicated. Black circles correspond to objects 
in the Serpens Cloud, while gray circles are background sources (see text for details). Note 
the different scales in the two plots. 



the theoretical median age of a molecular cloud (IBallesteros-Paredes fc Hartmannll2007l ). we 
cannot rule out from the current age distribution that several star-forming events have taken 
place in the observed Serpens cloud. This is also consistent with the studied region in Serpens 
being larger than other c oeval star-forming clu sters, where less di sperse age distr ibutions have 
been found, e.g. IC 348 (iMuench et al.l 120071 ) or Chamaeleon I (iLuhmanl 120071 ). 



Some basic tests for environmental trends were performed using our current dataset 
consisting of 78 stars. No strong correlations between the main parameters derived from the 
stellar spectra (ages and masses) and their positions in the cloud were found. This analysis 
will be revisited once the entire sample of 150 YSOs will be available. 

Recent analysis suggests that the distance to Serpens may be as low as 193 ± 13 pc 
(Jens Knude, private communication), rather than the distance of 259 pc assumed here. 
The smaller distance would imply a notable increase in individual ages due to the dec rease 
in luminosi t ies, w ith median ages of 15.9 and 10.1 Myr with the iBaraffe et al.l (1l998l ) and 
Siess et al.l (l2000l ) tracks, respectively. Median masses are not so strongly affected due to 
the mass tracks being almos t vertical in the low- mass regime, yie l ding m edian masses of 0.8 
and 0.6 Mq according to the lBaraffe et al.l (119981 ) and lSiess et al.l (120001 ) tracks, respectively. 
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Fig. 8. — His tograms of masses a nd ages for the sample of YS Os in Serpens, derived from 
the models of iBaraffe et al.l (Il998l ) (top) and lSiess et al.l (120001 ) (bottom). 



7. Accretion Based on Ha Emission 



A relation between th e Ha line-to-continuum ratio and evolution of its emi t ting source has 



been widely explored (IBertout 



1989 



Hartmann 1998, MuzeroUe et al. 2003. White fc Basri 



20031 . Istahler fc Pallalbooi iNatta et al.ll2004f ). It is argued that the strength of the Ha line 



decreases with evolutionary stage. A similar decline is seen in the strength of the prominent 
Ca II H and K lines (3968 and 3934 A, respectively) and in the surface rotational velocity, 
re-enforcing the motivation for a chromospheric origin for main-sequence Ha emission. How- 
ever, ac tive T Tauri stars have Ha fluxes much larger than those predicted by chromospheric 
models (jStahler fc Pallall2004l ). meaning that for these objects Ha is unlikely to originate pri- 
marily in a chromosphere. Mass accretion from a circumstella r disk is thought to be respon - 
sible for the excess emission of classical T Tauri stars (CTTS) (ILynden-Bell fc Pringldll974j ). 
The inner disk is disrupted by the stellar magnetic field, resulting in magnetospher ic accre- 
tion, where disk material is channelled along the magnetic field lines into the star (Koenig] 
19911 ). Prominent emission lines originate in the free-falling gas (jCalvet fc Hartmannlll992l ). 
Ha is one of the strongest emission lines in CTTSs, and has been extensively studied as an 
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indicator of accretion (e.g. Iwhite fc Basrilbool MuzeroUe et al.ll2003l . iNatta et al.ll2004h . 



We use the Ha line, covered by our optical spectra and present in emission in 35 objects 
of this sample, to obtain an indication of accretion activities according to two different 
methods from the literature: the Ha equivalent width (EW) and its full width at 10% of the 
peak intensity. The non-detection of Ha is not uncommon for pre-main-sequence stars. The 
sources could either be more evolved WTTSs (see § 17.11) . or th e non-detection could also be 
explained by the known temporal variability of T Tauri stars (jStahler fc Fallal l2004j ). 



7.1. Ha Equivalent Width 



According to lWhite fc Basril (120031 ) . there is no unique EW[Ha] value to distinguish between 
accreting CTTSs and weak-line T Tauri stars (WTTS). For instance. Ha emission from 
equally saturated chromospheres will appear much more prominently in a late-M star than 
in an early-K star due to the former's substantially diminished photospheric continuum near 
6500 A. 



Taking into account this spectral type dependence, IWhite fc Basril (120031 ) proposed that 
a T Tauri star is classical if EW[Ha] > 3 A for K0-K5 stars, EW[Ha] > 10 A for K7-M2.5 
stars, EW[Ha] > 20 A for M3-M5.5 stars, and EW[Ha] > 40 A for M6-M7.5 stars. These 
values were determined empirically from the maximum values of EW[Ha] for non- veiled 
T Tauri stars within each spectral type, from high-resolution spectra. Stars with values of 
EW[Ha] below these levels are not necessarily WTTSs, for which confirmation depends upon 
the Li abundance or other youth indicators, such as variability. 

We have determined the equivalent width of Ha with the usual fitting of Gaussian 
profiles to the lines and the results are given in Table |H The Ha feature lies near the 6567 A 
TiO band head, a strong feature in mid- and late-M stars. In low resolution spectra, the edge 
of this band becomes blended with the Ha emission feature, leading to an underestimate 
of the continuum on the redward side and thus an overestimate of the EW[Ha]. Also, self- 
absorption is missed in these unresolved data. For these reasons, our EW[Ha] values are 
likely somewhat overestimated. Nonetheless, we identify 35 stars with large EW[Ha] to be 
classified as CTTSs. The remaining four stars present marginal detections. 



7.2. Full Width of Ha at 10% 



Due to the uncertainties in the determination of EW[Ha] in low resolution spectra, IWhite fc Basri 
(120031 ) proposed the full width of Ha at 10% of the peak intensity (Ha [10%]) as an alterna- 
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tive to distinguish CTTSs from WTTSs. They found that stars with Ha [10%] greater than 
270 km s~^ are CTTSs, independent of spectral types. 

We fitted Gaussian profiles to the Ha lines (shown in Figure [H]) in order to determine 
the Ha [10%]. The errors were calculated by propagating the error on the spectral resolution 
and the error of the gaussian fit to the Ha line. The measured FWHM of each profile was 
deconvolved assuming a gaussian instrumental profile. Because of the low resolution of our 
spectra (183.6 km s'^ for WYFFOS, 126 km s'^ for DOLORES, and 100 km s"^ for CAFOS), 
very narrow profiles (e.g. objects 106 and 150) are unresolved. 

In spite of the low resolution, the Ha profiles found are broad enough in most cases to 
allow the extraction of velocity information, although with low accuracy. These results are 
shown in TableHl According to this criterion, 33 objects are CTTSs, while the remaining four 
stars are consistent with being CTTSs, given their uncertainties, showing a good agreement 
between both methods described. 




-1000 -500 500 1000 -500 500 1000 



AV (km/s) 

Fig. 9. — Continuum-subtracted profiles of Ha emission lines. The solid red lines are the 
Gaussian fit to each profile, dashed red lines mark the center of the Gaussian. 
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Studying accretion properties in low mass objects, iNatta et al.l (120041 ) found that Ha can be 
used not only as a qualitative indicator of the accreting nature of low mass objects, but also 
to obtain a quantitative estimate of the mass accretion rate Mac- They found an empirical 
relationship between the mass accretion rate, Mac (from Ha profile model fittings or veiling 
measurements) and Ha [10%] from high resolution spectra (Figure 3 in their paper): 



logMac = -12.89(±0.3) + 9.7(±0.7) x 10"3Ha[10%] (1) 

where Ha[10%] is in km s~^ and Mac is in M0yr^^ We use Eq. [1] to derive Mac for the 
objects showing Ha in emission (see Table Hj). 

The Mac values derived here (presented in Figure [TOl) are largely consistent with the 
range typical for T Tauri stars. Exception are objects 109, 50, 7 and 68 (Ha [10%] = 175, 
762, 805 and 966 km s~^, respectively). These Ha [10%] values are out of the range over which 
relationship ([T]) was calibrated, which may result in erroneous accretion rates (shaded areas 
in Figure fTOl) . The distribution of mass accretion rates versus ages shows a large scatter, 
although for the 7 strongest accreters (logMac > —7) a weak trend is seen, in the sense of a 
decreasing accretion rate function of stellar age. 




oLL . , . . , , . . , L_ 

-10 -8. -6 

log (Mo yr-') 



Fig. 10. — Histogram of mass accretion rates derived from Eq. ([T]) for the sample of sources 
in Serpens with Ha in emission. The shaded areas show obje cts out of th e rang e for which 
the relationship between Ha [10%] and Mac was calibrated by iNatta et al.l (120041 ). 



8. Conclusions 



We have presented a spectroscopic survey at optical wavelengths designed to determine 
spectral types and confirm the pre-main sequence nature of a sample of young stellar ob- 
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ject candidates i n the Serpens Molecu lar Cloud, selected on the basis of the "c2d" Spitzer 
Legacy Program (IHarvey et al.ll2007bl ). This sample will subsequently be used to study disk 
evolution, making use of Spitzer infrared photometry and spectra, obtained by our group. 



We have determined spectral types for 78 stars in the cloud, finding mostly K and M 
stars (88%), a late- type young stellar population. However, the high optical extinction 
toward this region makes it impossible to detect about half of the original sample of 
150 stars with optical spectrosco py. Near-IR spectroscop y will be used to classify the 



remaining objects in the sample (jOliveira et al. 



m prep 



Comparison of extinction values determined from our optical spectra with those cal- 
culated from the contribution of the cloud at each point (averaged over a 5' radius) is 
a powerful method for identifying candidate background sources with infrared excess 
from Spitzer selected samples of YSO candidates. The effectiveness of this method has 
been confirmed for two objects by follow-up IRS spectra. 

Effective temperatures and luminosities were derived and the objects were placed in a 
H-R diagrara , comp ared with theoretical i sochrones and mass tracks from models of 
Baraffe et al.l (119981 ) and ISiess et al.l (120001 ). 20 objects are too luminous to be at the 
distance of Serpens. These objects have a nearly one-to-one correspondence with those 
identified by the extinction method and are therefore very likely background objects. 
They amount to 25% of the detected sample, a relatively large amount likely due to 
the position of the Serpens Molecular Cloud at very low Galactic latitude. Two of 
these background sources were confirmed to be post-AGB stars through IRS spectra. 



The theoretical models by iBaraffe et al.l (119981 ) and ISiess et al.l ( l2000l ) imply a popu 
lation of YSOs concentrated in the age range between 1 and 15 Myr, st rongly peaked 
at 2 - 6 My r. The med i an ag e was found to be 7.5 and 4.7 Myr with the iBaraffe et al. 



( 119981 ) and ISiess et al.l (120001 ) tracks, respectively. Individu al masses raii g e from 0.2 
to 2.2 M(?t, with median values of 0.8 and 0.6 Mq for the IBaraffe et al.l (119981 ) and 
Siess et al.l (120001 ) models, respectively. 



The optical spectra cover the Ha line, an important indicator of accretion. We explored 
this relationship through two different approaches: qualitatively, the accreting nature 
of low-mass objects has been determined by either the equivalent width of Ha, or its 
full width at 10% of Ha peak intensity. This confirms 37 objects (or 55% of the YSO 
sample) to be actively accreting objects classified as new classical T Tauri stars in 
Serpens. The quantitative estimate of Mac based on the full width of Ha at 10% of the 
peak intensity yields a median of ~ 10~^ M0yr~^, with a broad distribution of values. 
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Table 1: Observation log 



Date 


Telescope + 


Wavelength 


Resolution 


Spectral 


Object ID (as in Table [2]) 




Instrument 


Blaze (A) 


(A) 


Coverage (A) 




2006 May 4 


WHT+WYFFOS 


7000 


4.0 


3000 


16, 34, 44, 47, 48, 87, 89, 91, 93, 96, 97, 
113, 119, 123, 125, 128, 130, 145, 148 


2006 June 29, 30 


TNG+DOLORES 


7000 


2.8 


1600 


2, 3, 4, 6, 35, 52, 70, 76, 85, 86, 108, 
118, 120, 122, 131, 139, 146, 150 


2007 June 9, 10, 11 


WHT+WYFFOS 


7000 


4.0 


3000 


2, 4, 5, 7, 14, 16, 18, 22, 23, 29, 30, 35, 
36, 40, 43, 52, 53, 55, 60, 61, 62, 66, 70, 
71, 76, 77, 81, 82, 84, 85, 86, 87, 88, 92, 
93, 96, 98, 99, 106, 110, 112, 113, 114, 
117, 118, 119, 120, 121, 122, 124, 125, 
127, 130, 131, 138, 139, 141, 142, 146, 
148, 149, 150 


2008 June 29, 30 


Calar Alto+CAFOS 


7750 


2.2 


3500 


1, 41, 58, 70, 82, 93, 99, 115, 120, 130, 
139, 141, 149, 150 
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Table 2. Stellar parameters for the YSO candidates in Serpens 



# 




c2d ID 




Spectral Type 












SSTc2dJ 




(isubtype) 


(mag) 


(mag) 


(mag) 


1 


1 


18275383-0002335 


-1.26 


K2±3 


1.0 


2.98 


9.0 


2 


2 


18280503+0006591 


-2.29 


M4±l 


6.5 


1.98 


2.0 


3 


3 


18280845-0001064 


-0.70 


M0±1 


0.0 


2.85 


2.0 


4 


4 


18281100-0001395 


-2.13 


M5±l 


8.0 


2.85 


6.0 


5 




18281315+0003128 


-1.80 


M3±3 


14.0 


4.83 


9.0 


6 


5 


18281350-0002491 


-1.08 


K5±2 


2.0 


5.33 


3.0 


7 


6 


18281501-0002588 


-0.05 


M0±1 


2.0 


5.33 


7.0 


14 


13 


18282143+0010411 


-1.44 


M2±2 


3.5 


2.66 


3.0 


16 


15 


18282432+0034545 


-2.01 


M0±1 


8.0 


3.41 


6.0 


18 


16 


18282738-0011499 


-1.98 


M5±4 


10.0 


4.78 


11.0 


22 




18283000+0020147 


-1.60 


M4±3 


14.0 


4.74 


12.0 


23 




18283736+0019276 


-2.19 


M7.5+1 


0.0 


7.59 


11.0 


29 


25 


18284481+0048085 


-1.01 


M2±3 


3.0 


4.29 


6.0 


30 


27 


18284497+0045239 


-1.35 


Ml±l 


2.5 


7.55 


2.0 


34 




18284828-0005300 


-2.43 


M5.5+1.5 


11.0 


4.92 


10.0 


35^ 




18284938-0006046 


-2.01 


M5±l 


12.0 


5.72 


13.0 


36 


32 


18285020+0009497 


-0.24 


K5±4 


5.0 


5.40 


10.0 


40 


36 


18285249+0020260 


-0.15 


M7±4 


8.5 


10.00 


12.0 


41 


37 


18285276+0028466 


0.06 


K2±3 


0.0 


19.20 


10.0 


43 


39 


18285395+0045530 


-1.16 


MO.5+3 


3.5 


8.49 


7.0 


44 


40 


18285404+0029299 


1.35 


M6±5 


0.0 


15.60 




47 


43 


18285489+0018326 


0.88 


M5±4 


13.5 


14.20 




48 


44 


18285529+0020522 


-0.29 


M5.5+3 


3.0 


6.86 


14.0 


52 


48 


18285808+0017244 


-2.11 


G3±5 


3.5 


9.87 


6.0 


53 


49 


18285860+0048594 


-1.06 


M2.5+3 


2.5 


3.37 


6.0 


55 


51 


18290025+0016580 


-1.01 


K2±3 


4.5 


7.88 


12.0 


58 


54 


18290088+0029315 


-0.42 


K7±2 


3.5 


19.00 


5.0 


60 


56 


18290122+0029330 


-0.39 


MO.5+0.5 


3.0 


19.00 


7.0 


61 


58 


18290175+0029465 


-0.84 


M0±1 


2.0 


19.00 


5.0 


62 


59 


18290184+0029546 


-0.77 


K0±7 


4.0 


19.00 


8.0 


66 


62 


18290393+0020217 


-0.82 


K5±2 


1.5 


3.53 


7.0 


70 


65 


18290575+0022325 


-2.12 


A3±3 


2.5 


4.17 


5.0 


71 


66 


18290615+0019444 


-1.52 


M3±l 


2.0 


4.06 


6.0 


76 


72 


18290775+0054037 


-1.06 


Ml±l 


2.0 


2.32 


6.0 


77 


73 


18290808-0007371 


-2.22 


M4±l 


8.5 


3.10 


8.0 


81 


77 


18290980+0034459 


-0.97 


M5±4 


1.0 


11.00 


15.0 
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Table 2 — Continued 



# 


ID"* 


c2d ID 


Ct2/:tm-24(Lnn'' 


Spectral Type 


AyOhs" 


AyCloud'^ 


AyFina 






SSTc2dJ 




(isubtype) 


(mag) 


(mag) 


(mag) 


82 


78 


18291148+0020387 


-1.34 


MO+2 


2.0 


3.30 


6.0 


84 


81 


18291407+0002589 


-2.19 


M3+3 


9.0 


4.74 


9.0 


85^= 




18291477-0004237 


-2.29 


M3+1 


8.0 


2.71 


9.0 


86 


83 


18291508+0052124 


-1.89 


M5.5+1.5 


8.0 


2.51 


8.0 


87 


84 


18291513+0039378 


-1.50 


M4+3 


2.5 


7.95 


7.0 


88 


85 


18291539-0012519 


-1.72 


MO. 5+2 


7.0 


2.51 


7.0 


89 


86 


18291557+0039119 


-0.56 


K5+5 


2.0 


10.20 


12.0 


91 


88 


18291617+0018227 


0.45 


K7+3 


1.0 


4.97 


18.0 


92 


89 


18291969+0018031 


-1.23 


MO+1 


2.5 


3.51 


7.0 


93 


90 


18292001+0024497 


-2.04 


M2+3 


8.5 


3.17 


9.0 


96 


94 


18292184+0019386 


-0.99 


Ml+1 


2.0 


3.16 


7.0 


97 




18292250+0034118 


-1.65 


M2+2 


3.0 


10.10 


10.1 


98 




18292253+0034176 


-2.48 


A3+5 


5.5 


10.10 


5.0 


99 


95 


18292616+0020518 


-2.26 


M4+2 


13.0 


3.34 


8.0 


106 


102 


18292927+0018000 


-1.24 


M3+1 


2.5 


3.48 


7.5 


108 


105 


18293254-0013233 


-2.58 


K5+5 


0.0 


2.51 


1.0 


110 


107 


18293319+0012122 


-1.81 


M6+2 


8.5 


3.90 


5.0 


112 


109 


18293381+0053118 


-2.06 


M7+4 


1.5 


4.59 


8.0 


113 


110 


18293561+0035038 


-1.53 


K7+1 


2.0 


12.60 


4.0 


114 


111 


18293619+0042167 


-0.92 


F9+5 


7.0 


6.71 


11.0 


115 


114 


18293672+0047579 


-1.08 


MO.5+2 


6.5 


5.03 


8.0 


117 


118 


18294020+0015131 


0.87 


K2+5 


0.5 


3.40 


2.0 


118 




18294067-0007033 


-2.53 


M4+1.5 


8.0 


1.41 


8.0 


119 


119 


18294121+0049020 


-1.39 


K7+2 


4.0 


5.20 


5.5 


120 


123 


18294168+0044270 


-1.42 


A2+2 


5.0 


6.13 


8.5 


121 




18294301-0016083 


-2.18 


K7+5 


5.5 


3.33 


7.0 


122 


126 


18294410+0033561 


-1.68 


MO+1.5 


1.5 


7.35 


5.0 


123 


129 


18294503+0035266 


-1.26 


MO+1 


1.5 


7.02 


10.0 


124 


132 


18294725+0039556 


-1.32 


MO+1.5 


2.0 


10.40 


5.5 


125 


133 


18294726+0032230 


-1.31 


MO+1 


3.0 


6.59 


7.0 


127 


143 


18295001+0051015 


-1.68 


M2+1 


4.0 


4.18 


6.5 


130 


145 


18295041+0043437 


-1.32 


K6+1 


4.5 


8.69 


8.5 


131 


148 


18295130+0027479 


-2.42 


A3+3 


2.0 


3.34 


5.5 


138 


161 


18295322+0033129 


-1.99 


MO+3 


10.0 


6.86 


11.0 


139 


165 


18295422+0045076 


-2.35 


A4+4 


4.0 


6.95 


6.5 


141 


169 


18295531+0049393 


-0.74 


A3+3 


3.5 


3.02 


9.0 
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Table 2 — Continued 



# ID- 


c2d ID 


a2Mm-24^un'' Spectral Type 


AvOhs'^ 


AvCloud'^ AyFinal 




SSTc2dJ 


(isubtype) 


(mag) 


(mag) (mag) 



142 


172 


18295592+0040150 


-0.85 


M4±2 


3.5 


5.71 


3.5 


145 


178 


18295714+0033185 


-0.82 


G2. 5+2.5 


3.5 


5.19 


10.0 


146 


182 


18295772+0114057 


0.44 


M4+2 


4.0 


14.05 


3.0 


148 


206 


18300178+0032162 


-1.25 


K7+1 


2.0 


5.00 


7.0 


149 


210 


18300350+0023450 


-1.22 


MO+1 


0.0 


2.25 


3.5 


150 


222 


18300861+0058466 


-0.70 


K5+3 


0.0 


2.79 


5.0 



^From lHarvev et al.l |2007b[) 



'^Obtained from a li near fit to the logarit hm of the fluxes between the 2MASS K (2/im) and the 
MIPSl (24/im) bands |Harvev et al.ll2007bl ) 

'^Extinction found in this work 

'^Extinction found by "c2d" over 5' region 

''post-AGB stars confirmed with Spitzer IRS spectra 
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Table 3. Stellar parameters for the YSO candidates in Serpens 



# 


Teff 




Ages 


MassB 


Ages 


Masss 




(K) 


(Lo) 


(Myr) 


(Mo) 


(Myr) 


(Mo) 



1 


4900 +11° 


2.01+2-39 


11.60+13-91 


1, 


.561°; 


63 


6.10+iii°° 

— 4. 10 


1, 


,57+°- 

— U. 


,43 

44 


2t 


3270 ±145 


5.24+6-26 

— Z.i)0 


















3 


3850 


0*51_Q 29 


5.42ti°6f 


0, 


.921°; 


13 
.16 




0, 


•571^; 


,22 
.11 


4t 


3125 


1.30+n-7^ 


















5t 


3415 ttll 


14.78li^3f 


















6 


4350 +fcPn 


1 13+^-35 

— 0.64 


5.98+13219 


1, 


,381°; 


18 

.04 


3.41+8-^^ 

— Z. iU 


1, 


,13+° 


,19 

.0 i 


7 


3850 


0.23+0-" 

— U.io 


16.25+?^-98 


0, 


.811°; 


,11 

.15 


8.26+^1,^ 

— 0.4o 


0, 


,59+° 

— U. 


,19 
, 12 


14 


3560 ±290 


0.19t0-?2 


6.87+1^9™ 


0.52±0 


.25 


4.37+^-90 


0, 


,37+°- 


.22 

12 


16 


3850 +11° 




1.56+3-^^ 

— I.OD 


1, 


.03+°- 

— U. 


16 

1 7 


0-52l2;i9 


0, 


,55+°- 

— U. 


^20 
no 


18t 


3125 

— 1 ZD 


25.26+?"-j^^ 


















22t 


3270 +tll 


81.97+^r?? 


















23t 


2795 


1 Q2+2-29 

— 1 .Uo 


















29 


3560 +111 

41:00 


n+n-14 

U . -L J- Q 


-|9 r,i +62.38 


0, 


■5oii;' 


16 

.00 


7.80+f8|4 

— 4.00 


0, 


,36+E} 


,29 

1 Q 


30 


3705 ±145 


0.30+°-36 

— U.i i 


6.57+13395 

— 4.00 


0, 


.711°; 


15 
.16 


3.41+?i^ 

— i.oo 


0, 


,461°; 


,12 

.09 


34t 


3057 +177 


3-24+?-S 

— L.oZ 


















35''' 


t 3125 tlfn 


306.47tf7^-Ig 


















36 


4350 +llr 

— d45 


1-09+n-^? 

— U.bl 


6.34+19,20 

" " — 5.33 


1, 


.361°; 


,20 
.51 


3.68+i*4f 

— o.4d 


1, 


.131E!: 


,03 
.67 


40 


2880 +fl 


0.23+°-2^ 

— U.io 










2.86+2-86 

— z.ou 


0, 


.311^: 


,31 
.30 


41 


4900 +ll°o 


0.01±0.01 


















43 


3777 +™ 


0.33+°-?^ 

—0.18 


7.47+19-f 
' ' —6.17 


0, 


.801°; 


11 

.44 


3.79+"n^^ 
" • ^—2.02 


0, 


,521°; 


.39 
.23 


44b 


2990 till 




















47b 


3125 +^25 




















48 


3057 


U.OD_o 20 










2.4712;^^ 


0, 


,351°; 


.35 
.30 


52 


5830 


'"-2.68 


8.70ti^„f 


2, 


.53+?- 


,64 
,16 


12.00l9;697 


1, 


,501- 


,46 
,28 


53 


3487 1^37 




'-'•.'■'^-0.07 


7 c 9+49.36 
' -^^-6.95 


0, 


.441°; 


,27 
,31 


r or+34.68 

^■•^"^-3.54 


0, 


,33lS; 


,35 
.17 


55 


4900 




9.45ti23f 


1, 


.601°; 


,53 
,45 


4.98l9;i3 


1. 


,661°; 


,41 
.53 


58 


4060 tlf. 


1 00+1-1^ 

l.UU_Q 5Q 


1 c+7.30 
"^•-■^^-2.20 


1, 


.181°; 


15 

,34 


2.39l3;93^ 


0, 


,761°; 


,37 
.30 


60 


3777 ±72 




r 70 + 11.95 

'"^-3.81 


0, 


,821°; 


,08 
,09 


2.89l?;°7l 


0, 


,521°; 


,06 
.05 


61 


3850 tl^^ 


-1 Q7+1-63 

-0.77 


i.24t?-;i 


1, 


,051°; 


17 

,16 


2.07±0.41 


0.66±0 


.09 


62 


<^9^C\ +580 


6.70t^;f7 


2.08t^;i? 


2, 


.191J; 


17 
,56 


07+4.94 

-3.63 


2. 09+? 


,49 

.52 


66 


4350 tll'o 


1.87t?;23 


A.2Q+_r,i 


1. 


.441°; 


,25 
,19 


9 1 q+3.49 

Z.iO_2.i3 


I.I2I- 


,42 
,39 


70 


8720 t^™ 


12.64+1^1^ 










8.18±2.21 


2.09±0 


.19 


71 


3415 ±145 


o.20t°;?t 


2.86i?;90 


0, 


,381°; 


14 
.11 


9 04+3-96 
^•^^-0.93 


0.31±0 


.06 


76 


3705 ±145 


0.21^0;?^ 


9.62121^62 


0, 


.681°; 


11 

.16 


c 40+11.80 


0.46+° 


.13 
.09 


77t 


3270 ±145 


9 00+3.44 


















8lt 


3125 +_lll 


22.08126-35 
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Table 3 — Continued 





Toir 


L* 


Ages 


Mass B 


Ages 


Masss 




(K) 


(Lo) 


(Myr) 


(Mo) 


(Myr) 


(Mo) 



82 


3850 


+355 
-290 


0. 


iq+0.15 
— 0.07 




0, 


•721°: 


.13 
.22 


19.65lf^.^« 


0.58ir2l 


84t 


3415 


+435 
-425 


8. 


qq+9^94 

■ OO ^ gg 














85a,t 


3415 ±145 


33 


00+40.37 
•0^ — 19.03 














86t 


3057 


+212 
-177 


22 


r9+26.'89 
•3^-12.67 














87 


3270 


+435 
-390 


0. 


90+0. 28 

■ ^0 Q "1^2 




0, 


.28+;5' 


,41 
21 


C-7+2.13 

Zj .0 1 2 57 


0.25 Q Qg 


88 


3777 


+355 
-290 


1. 


4n+l>7 

■*°-0 83 




0, 


.97l"o' 


:3i 

14 


^- '"-0.27 




89 


4350 


+900 
-790 


0. 


4t;+0.54 
■^"-'-0.25 


21.77tt9'5f 


1, 


,oolEJ 


,32 
42 


13.98l3|23 




91 


4060 


+530 
-500 


2. 


00+2.79 
'"^'^— 1.31 


0.8212-26 


1, 


,28ll]' 


59 
24 


2.53±0.20 


1.13±0.37 


92 


3850 


+210 
-145 


0. 


r,n+0.34 

"-0 16 


10.821^3^60 


0, 


.851°: 


:io 
,15 


r 09 + 14.61 
J.o.i_3 42 


0.58l°;22 


93t 


3560 


+500 
-435 


4. 


oq+5.24 
'^''-2.47 














96 


3705 ±145 


0. 


91 +0.26 
■^'^—0 12 


9.461^1^21 


0, 


.681° 


,11 

16 


r 90 + 11.50 

j.z;o_2 7g 


0.46l°:J^ 


97 


3560 ±290 


0. 


-if;+0.19 
■^'J-O 09 


8.011^^8^ 


0, 


.5itg: 


,24 
,25 


4.94l^;f7 


ri 07+0.22 
U.O 1 _Q 12 


98 


8720 


+3180 
-1140 


5. 


44+6-49 
■^^-3.06 










3.7911:^9 


9 70+2.73 
^. 1 o_2 73 


99t 


3270 


+290 
-280 


17 


97+20.61 
' —9.72 














106 


3415 ±145 


0. 


1 7+0.20 
■J- ' _0 09 




0, 


.371°: 


,14 
,10 


3.631-^:22 


0.30±0.06 


108t 


4350 


+900 
-790 


14 


^0 + 17. 28 
•4°-8 15 














110^ 


2990 ±280 


4. 


87+5.81 
' -2.74 














112t 


2880 


+535 
-880 


2. 


79+3.24 
' ' —1 53 














113 


4060 


+ 145 
-210 


0. 


C9+O.6I 
■"•^—O 29 


8.641^^8^° 


1, 


.061°: 


,13 
,22 


4.78l^:'6°i 


0.79l°:^6 


114 


6115 


+475 
-315 


1. 


nq+2.25 
■"•' — 1 06 










19.37±3.64 


1.30±0.09 


115 


3777 


+355 
-290 


0. 


9^+0.29 
■^^-0.14 


10.28t7^2f 


0, 


.771°: 


,08 
,32 


5.701^:^2^ 


'-'•"J^-0.17 


117 


4900 


+736 
-840 


0. 


09+0.02 
■"^-0.01 














118t 


3270 


+217 
-212 


19 


f.7+23.47 
-11.06 














119 


4060 


+290 
-355 


0. 


1 c+0.18 
■-^"J-0.08 




0, 


,731°: 


,20 
,07 


36.46126^6 


70+° " 

^■'"-0.25 


120 


8970 


+550 
-510 


10 


^0 + 12.51 
•40-5.89 










7 00 + 7.22 
/ .ZiZi y -^-j 


9 iq+2.19 
z.iy_2.i4 


121 


4060 


+840 
-790 


0. 


qQ+0.47 
■"^^-0.22 




0, 


.991- 


15 
,69 


7 00+40.69 
( .00 5 Y3 


P) 0-1+0.06 
L'-«^-0.55 


122 


3850 


+282 
-217 


0. 


54+0.64 
■^^-0.30 




0, 


■931°: 


,17 
,25 


2-681^:°^, 


0.571°:?^ 


123 


3850 


+210 

-145 


0. 


qc+0.42 
■'-'^-0.20 


8.54li^^« 


0, 


,881°: 


,11 
,15 


4 04+10.41 

^•'J+-2.43 


0.581°:?? 


124 


3850 


+282 
-217 


0. 


17+0.21 
■-^'-0.10 


23.46+g-^^ 


0, 


,761°: 


,12 
,17 


12.65+f7f 


0.591°:!^ 


125 


3850 


+210 

-145 


0. 


90+0.34 

■'^°-0.16 


10.98tl%f 


0, 


,851°: 


,10 
,15 


5.89l^if 


0.581°:?^ 


127 


3560 ±145 


0, 


04+0.41 
■"J+-0.19 


9 07+5.76 
^■^'-1.83 


0, 


,56i;j: 


,17 
,13 


2.62l2-^^ 


0.381°:°? 


130 


4205 


+385 
-355 


0. 


70+0.87 
■ '"^-0.41 


7.88116^^1 


1, 


,191°: 


,18 
,27 


4 30+12-07 

^•"^^-2.68 


n Qr,+0.20 
*J-y^-0.38 


131 


8720 


+800 
-695 


13 


7C + 16.42 
• ' ^-7.74 










"^■^"-6.00 


1 04+0-42 
-■-•^+-0.18 


138t 


3850 


+500 
-435 


3. 


70+4.46 
■ '^-2.10 














139 


8460 


+ 1060 

-880 


17 


qfi+21-44 
•^"-10. 10 










QQ+40.03 
^•^^-8.80 


9 (.1+0.51 
^■^^-0.27 


141 


8720 


+800 
-695 


38 


on+46.31 
.OU_2i.83 










4 51 +5-16 
^•^■'■-1.84 


9 CI +0.54 
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Table 3 — Continued 



# 


Teff 




Ages 


Masss 


Ages 


Massg 




(K) 


(Lo) 


(Myr) 


(Mo) 


(Myr) 


(Mo) 



142 


3270 


+290 
-280 


-0.04 


5.701^9^ 


0.24 


145 


5845 


+ 185 
-75 


3.15 ^ 


180.47l^^^-J? 


2.01 


146 


3270 


+290 
-280 


n 97+0.33 
-0.15 


l-"O-1.05 


0.28 


148 


4060 


+ 145 
-210 


^■^^-0.29 


o 54+18.89 
°-"-'^-5.82 


1.06 


149 


3850 


+210 
-145 


^■^^-0.08 


9Q cq+66.82 


0.74 


150 


4350 


+550 
-500 


1 iq+1.35 
-'^■^'-'-0.63 




1.37 



0-26 R n 91 +013 

Q 0-UO_2.96 'J-^^-O.lO 

+ 1-34 1fi/1S+44.83 1 qq+0.36 

-0.81 l"-*°-7.78 1"J«J-0.17 

+0.25 9 qr + 1.75 n 9r+0.12 

0.13 ^■'-'^-59.44 '-'■^^-0.04 

+0.13 A 79+13.78 r> 7Q+0.16 

-0.22 ^■''^-2.62 "-'^-0.22 

+0.13 17QR+48.16 n c^Q+n io 

-0.17 -'^'■'-'"-12.23 '-'•'J^-0.13 

+0.18 q /ic+9.21 1 -1 q+0.12 

-0.35 '-'■^"-2.69 ^■-'^'-'-0.57 



'^post-AGB stars, confirmed with Spitzer IRS spectra 
''Objects without 2MASS J, H and K magnitudes 
^background objects 



Note. — Subsript B corresponds to values derive d from the 
models, while S corresponds to values derived from lSiess et al. 



Baraffe et all |l998l ) 
(l2000h 
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Table 4. Ho; equivalent widths and line widths, and mass accretion rates for YSOs in 

Serpens 



# 


(6563 A) 


CTTS 
(EW[Ha]) 


Ha[10%] 
(km/s) 


CTTS 
(Ha[10%]) 


log Mac 

(Mo/yr) 


i 


4.0 


/ 
V 


Add J_ OO 


/ 
V 


-o.o4 ± u.yo 


Q 

o 


01. 


/ 
V 


ooy =t oy 


/ 
V 


7 AA J- n 7p; 
- ( .00 =t U. /O 


D 




/ 
V 


QHK 1 1 1 n 

oUo =t iiU 


/ 
V 


no -U 1 O/i 
-O.Uo it l.z4 


/ 


00.0 


/ 
V 


404 zL o 


/ 
V 


QO -U n or; 
-o.oy zt U.yo 


1 A 


OQ O 

zo.y 


/ 
V 


AAQ -U 110 
445 =t iiZ 


/ 
V 


K/l -U 1 17 
-O.04 zt i.i / 


1 

io 














iU4.y 


/ 
V 


QQQ J- 0/1 
ooo =t o4 


/ 
V 


o AA J- n on 

-y.oo zt u.yu 


oU 


1 1 o 

ii.z 


/ 
V 


OQA J- A1 

zo4 zL Oi 


/ 
V 


1 n 1 Q J- n AO 

-iu.io zt u.oy 


OD 




/ 
V 


Q/in -U 000 
o4U =t zyz 


/ 
V 


O F^O -U OA 

-y.oy zt z.oo 


An 
4U 












A 1 


-U. / 










A Q 

4o 


1 O 1 


/ 
V 


0/1 1 J- C\A 

z4i zL y4 




1 o c j_ n 07 

-iu.oo zt u.y/ 


A A 
44 












4 ( 












A Q 


OOZ.O 


/ 
V 


7A0 -U 1 07 
/ Oz zt iz / 


/ 
V 


-O.OU zt i.o / 


52 


-2.2 










53 


92.7 




529 ± 112 




-7.76 ± 1.19 


DO 


110 
ii.Z 


/ 
V 


/1 1 Q _i_ 7n 

4io =t iVJ 


/ 
V 


00 -U n on 
-o.oo zt U.oU 


Do 


/ i.O 


/ 
V 


c:n7 -U AQ 
oU ( =t Oo 


/ 
V 


7 07 -U n 77 
- i .\) ( zt U. M 


OU 


44. o 


/ 
V 


04o =t iU4 


/ 
V 


7 AH -U 1 1 
- ( .OU zt i. iz 


Di 


zo.o 


/ 
V 


Q1 1 _L KT 

oil =t or 


/ 
V 


O 07 J- O A7 

-y.o/ zt U.O/ 


oz 












DO 


1 07 A 


/ 
V 


OCA _|_ 1 OQ 

yoo =t loo 


/ 
V 


Q pio -U 1 p;q 
-o.Oz zt i.Oo 


70 
1 U 


-±0.1 










71 


68.1 


^/ 


490 ± 84 




-8.14 ± 0.93 


76 


-9.7 










82 


25.2 




656 ± 123 




-6.53 ± 1.31 


87 












88 
89 


22.2 


V 


532 ± 174 


V 


-7.73 ± 1.75 


91 


108.9 


V 


549 ± 65 


V 


-7.56 ± 0.80 


92 


24.8 


V 


478 ± 138 


V 


-8.25 ± 1.41 


96 


22.0 


V 


704 ± 168 


V 


-6.06 ± 1.73 


97 












98 


-7.8 










106 


15.8 


V 


175 ± 186 




-11.19 ± 1.83 



-35- 



Table 4 — Continued 



# 




CTTS 


Ha[10%] 


CTTS 


log Mac 




(6563 A) 


(EWpa]) 


(km/s) 


(Ha[10%]) 


(Mo/yr) 


113 


10.9 


V 


273 ± 73 


V 


-10.24 ± 0.79 


114 


- 










115 












117 


91.9 


V 


559 ± 152 


V 


-7.47 ± 1.55 


1 1 Q 


Si S 
O.O 




fi'^Q -1- OA 


/ 
V 


-D.oy =t i.UO 


120 


-11.0 










121 












122 


37.1 


V 


444 ± 76 


V 


-8.58 ± 0.85 


123 












124 


34.7 


V 


225 ± 64 




-10.71 ± 0.71 


125 


49.0 


V 


618 ± 95 


V 


-6.89 ± 1.06 


127 


5.6 




460 ± 213 


V 


-8.43 ± 2.11 


130 


123.2 


V 


590 ± 114 


V 


-7.17 ± 1.22 


131 


-9.6 










139 


-6.4 










141 


9.0 


V 


289 ± 100 


V 


-10.09 ± 1.03 


142 


116.6 


V 


367 ± 75 


V 


-9.33 ± 0.83 


145 


-1.5 










146 


23.4 


V 


318 ± 102 


V 


-9.80 ± 1.06 


148 


18.6 


V 


598 ± 122 


V 


-7.09 ± 1.29 


149 


226.1 


V 


502 ± 36 


V 


-8.02 ± 0.58 


150 


3.4 


V 


213 ± 198 




-10.82 ± 1.95 



'^Positive values denote emission 



